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I To describe large momentum distributions of charged hadrons observed 

fH , at RHIC, a diffusion equation in the three dimensional hyperbolic space is 

Q_i' introduced. 

5^' PACS numbers: 25.75.-9, 25.75.Dw, 02.50.Ey 

1. Introduction 

^ I We are much intersted in events with high multiplicity at RHIC, where 

' thousands of particles are produced per event. To describe such many parti- 

cle system, a sort of collective approach will be useful. One-particle rapidity 
or pseudo-rapidity distributions observed at RHIC are well described by the 
Ornstein-Uhlenbeck process with two or three sources pp. In this paper, we 
would extend our approach to include an expansion for transverse direction. 
We can derive the formula same with that proposed by Voloshin |2j from 
the three-dimensional Ornstein-Uhlenbeck process using the variables, lon- 
gitudinal rapidity and two transverse momenta. The fundamental solution 
is gaussian in these variables. Even if the formula is integrated over the 
azimuthal angle in the transverse momentum plane, the distribution is only 
slightly modified from the gaussian. 

However, observed pt distributions at RHIC have long tails compared 
with an exponential distribution, and cannot be described by the formula. 



* Presented at XXXIII International Symposium on Multiparticle dynamics, Septem- 
ber 5-11, 2003, Krakow, Poland 
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Therefore relativistic approach, in other words, a stochastic model in the 
three dimensional hyperbolic space is adopted for the production processes, 
and pt distributions observed at RHIC are analysed by the model. 



2. Diffusion equation in the three dimensional hyperbolic space 

For simplicity, we consider the diffusion equation in the geodesic po- 
lar coordinate system. In this system, 4 energy- momentum (-E',Pl,Pt) of 
particle with mass m is specified by radial rapidity p, polar angle 9 and 
azimuthal angle (p. Then energy and absolute value of momentum are ex- 
pressed respectively as 



E = m cosh p, P = Y pi + Pt^ = m sinh p. 
The metric in the geodesic polar coordinate system is written as 

The diffusion equation, 
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is taken as a model of particle production processes: After a collision of 
nuclei, particles are produced at the origin of rapidity space expressed by 
Eq.©. Then those particles propagate according to the diffusion equation 
((2). In the course of the space-time development, energy is supplied from 
the leading particle system (composed of collided nuclei) to the produced 
particle system. Number density of particles becomes lower and at some 
(critical) density particles become free. 

Approximate solution of Eq. ^ for Dt < < 1 is given as [H] 
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Equation Q) is the solution of the radial symmetric diffusion equation [^, 
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if C = (47rL't)~^/^e~^*. In the geodesic cylindrical coordinate system, 4 
energy-momentum is expressed by longitudinal rapidity y, transverse rapid- 
ity ^ and azimuthal angle (p. Energy is given as 

E = m cosh y cosh ^ = mx cosh y. (6) 
where mT = \Jp\ + m"^- From Eqs.(^ and (jH)), we have an identity, 

cosh p = cosh y cosh ^. (7) 

Radial rapidity coincides with transverse rapidity, p = ^, when y = 0. 

Therefore, we can analyse transverse momentum (rapidity) distributions 
at fixed longitudinal rapidity (y = 0,2.2,---), using Eqs.Q) and with 
parameters, C and a{t)'^ 

3. Analysis of pt distributions observed at RHIC 

Transverse momemtum distributions of charged hadrons, and h~ , 
observed by the STAR [5 and BRAHMS 6 collaborations are analysed. 
Data are taken at fixed pseudo-rapidity ij. Particles are not identified in 
these experiments. Therefore, we use mass m as a parameter in addition 
to C and cr'^{t). To estimate the values of parameters, only the statistical 
error is used. 

The results at ry = by the STAR collaboration are shown in Fig.l and 
Table 1. Estimated value of m decreases as the centrality cut increases, 
where as variance cr'^it) decreases. 

The results at = by the BRAHMS collaboration are shown in Fig. 2 
and Table 2. Estimated value of m increases as centrality cut increases. 
However the value is somewhat higher than that by the STAR collabora- 
tion at the same centrality cut. Variance decreases as the centrality 
cut increases. In this case the value is smaller than that by the STAR 
collaboration. 



centrality 


m 


C 


a{tY 


• /n.d.f 


00-05% 


0.57 


316.7 ± 5.2 


0.332± 0.001 


244.9/32 


05-10% 


0.56 


266.5 ± 4.7 


0.334± 0.001 


129.6/32 


10-20% 


0.53 


202.8 ± 3.1 


0.354± 0.001 


144.9/32 


20-30% 


0.52 


143.5 ± 2.5 


0.378± 0.001 


97.94/32 


30-40% 


0.50 


101.2 ± 2.5 


0.360± 0.001 


98.3/32 


40-60% 


0.44 


62.7 ± 1.1 


0.402± 0.001 


52.1/32 


60-80% 


0.37 


25.6 ± 0.5 


0.444± 0.001 


26.6/32 


min-bias 


0.32 


3.03 ±0.08 


0.457± 0.004 


16.1/29 



Table 1. Estimated parameters of pt distributions of — ^ — at r/ = in 
Au+Au collisions at ^snn=200 GeV 0. 
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lir" r 



BRAHNS (2006eV) 
Au+Au. (h*+hT/2 at ti=0 



00-1 OJ 
■=0 71 
1 0-201/5 
■=0 67 
20-40J/I0 
m-O 63 
40-601/15 
■=0 57 




P, (GeV/c) 



h++h- 



Fig. 2. pt distributions of 
= in Au-|-Au collisions i6| 





BRAHIS (200 GeV) 






Au+Au. h at 


T|=2.2 










00-101 
-■=0 80 
IO-20%/5 
■=0 84 
20-401/10 i 
■=0 81 
40-60«/l5 
■=0. 70 
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P, (GeV/o) 

Fig. 3. Pt distributions of h~ at rj 
2.2 in Au-f Au collisions |6j 



centrality 


m 


C 




Y'^ • /n.d.f 


00-10% 


0.71 


121.0 ± 0.7 


0.307± 0.001 


351.2/23 


10-20% 


0.67 


88.7 ± 0.6 


0.325± 0.001 


233.5/23 


20-40% 


0.63 


56.9 ± 0.7 


0.340± 0.001 


97.8/23 


40-60% 


0.57 


27.5 ± 0.5 


0.358± 0.002 


28.7/23 



Table 2. Estimated parameters of pt distributions of 



h++h- 



at = in 



Au+Au collisions at ^/si^7=200 GeV by the BRAHMS collaboration 0. 
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The results at ry = 2.2 by the BRAHMS collaboration are shown in Fig. 3 
and Table 3. Data are well described by Eq.(|ll). In this analysis, we assume 
y=2.2, and use Eq.(I7|). 



centrality 


m 


C 


a{tr 


Y'^ • /n.d.f 


00-10% 


0.80 


52286 ± 1575 


0.470± 0.002 


93.7/8 


10-20% 


0.84 


41669 ± 1561 


0.452± 0.002 


30.5/8 


20-40% 


0.81 


19364 ± 862 


0.469± 0.003 


68.9/8 


40-60% 


0.70 


6070 ± 418 


0.524± 0.005 


3.1/8 



Table 3. Estimated parameters of pt distributions of /i at ry = 2.2 in 
Au+Au colhsions at 200 GeV by the BRAHMS collaboration 0. 



4. Summary and discussions 

In order to analyse large pt distributions of charged hadrons observed 
at RHIC, a stochastic process in the three dimensional hyperbolic space is 
taken as a model of particle production process. The solution is gaussian- 
like in radial rapidity. 

Transverse momentum distributions at r/ = observed by the STAR 
collaboration and those at ?? = and 2.2 by the BRAHMS collaboration 
are well described by the formula, if the value of mass contained in rapidity 
variable is used as a parameter. 
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